Repulsive shield between polar molecules 
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We propose and analyze a technique that allows to suppress inelastic collisions and simultaneously 
enhance elastic interactions between cold polar molecules. The main idea is to cancel the leading 
dipole-dipole interaction with a suitable combination of static electric and microwave fields in such 
a way that the remaining van-der-Waals-type potential forms a three-dimensional repulsive shield. 
We analyze the elastic and inelastic scattering cross sections relevant for evaporative cooling of polar 
molecules and discuss the prospect for the creation of crystalline structures. 
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It is well known that an attractive van der Waals in- 
teraction dominates the microscopic interaction potential 
between ground state atoms and molecules at short dis- 
tances P| . As a consequence, ultracold atomic and molec- 
ular gases are prone to the formation of deeply bound 
molecular states via inelastic collisions. These inelastic 
collisions limit the lifetime of strongly interacting ultra- 
cold gases P, 0] ■ Controlling inelastic collisions is there- 
fore crucial for the creation of novel strongly correlated 
quantum degenerate systems, such as polar molecules. 
In this Letter, we demonstrate the possibility to engi- 
neer a three-dimensional repulsive interaction between 
polar molecules, which allows for the suppression of in- 
clastic collisions, while simultaneously enhancing clastic 
collisions. This technique may open up a way towards ef- 
ficient evaporative cooling and the creation of novel long- 
lived quantum degenerate gases of polar molecules. 

A special property of polar molecules prepared into 
the lowest rotational and vibrational state is a perma- 
nent electric dipole moment d, which gives rise to tun- 
able dipole-dipole interactions and offers a large po- 
tential for the creation of strongly correlated quantum 
phases 0, 0, 0, 01- Two routes are currently explored 
for the experimental realization of quantum degenerate 
polar molecules^ (V) trapping and cooling of preexist- 



ing molecules 0, S B EE ED, El and (ii) synthesiz 
ing polar molecules from a cold mixture of atomic gases 
EE El EE EE EE EE E! • While scattering properties of 
heteronuclear molecules with dipole-dipole interactions 
are currently theoretically explored 2^, 21, 22 1, it is ex- 
pected that inelastic collisions strongly increase for polar 
molecules compared to atomic systems due to the open- 
ing of additional decay channels. 

The main idea of our approach is to cancel the leading 
dipole-dipole interaction with a suitable combination of 
a static electric field and a continuous-wave microwave 
field: the former induces a dipole moment d z , while the 
latter drives an additional dipole moment d± rotating 



with frequency lj of the microwave field, see Fig.pjb). In 
analogy with magic-angle techniques in NMR [23|, the 
time-averaged interaction of the rotating dipole moment 
shows a negative sign allowing for a cancelation of the to- 
tal dipole-dipole interaction. The remaining interaction 
is tunable in strength with a repulsive van der Waals be- 
havior V et{ ~ (d 4 / HA) I Y 6 , where A is the detuning of 
the microwave field and r is the intermolecular separa- 
tion. The three-dimensional shield described here is thus 
purely repulsive. This is in contrast to the "blue shield" 
discussed in the context of atomic gases, which is attrac- 
tive at large distances P]. We find that the efficiency of 
the shield is determined by a single dimensionlcss param- 
eter 7 = d 2 m/h 2 r B with r B = (d 2 /B) 1 ' 3 , B the rotational 
energy, and m the mass of the molecule. Under optimal 
conditions, for large values of 7, inealstic collisions can 
be quenched for temperatures T < 0.01B. 
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FIG. 1: (a) Energy levels of H^l: the arrow indicates the 
microwave field, (b) Sketch for the cancelation of the dipole- 
dipole interaction: d z denotes the dipole moment induced by 
the static electric field, while the dipole moment d± is rotating 
in an orthogonal plane due to the microwave coupling. 



To describe the shield quantitatively, we consider two 
polar molecules prepared into the vibrational ground 
state. The interaction at large intermolecular distances is 
dominated by the dipole-dipole interaction, yielding the 
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where we have introduced the center of mass R = 
(ri + r 2 ) /2 and the relative coordinate r = ri — r 2 , with 
the corresponding momenta P and p, respectively, while 
r = |r| and r = r/r. Here, we are interested in po- 
lar molecules in the 1 £ electronic ground state, and the 
Hamiltonian for the internal structure takes the form [241 



HQ 



B J? 



dj • E dc — dj • Eac(i) 



(2) 



with the dipole operator d^ and the permanent dipole 
moment d. The first term describes a rigid rotor account- 
ing for rotational structure with the rotational energy B, 
while the last two terms describe the coupling to an ex- 
ternal static electric field E dc and microwave field E ac . 
While additional interactions with the nuclear spins are 
in general small and can be ignored, the analysis pre- 
sented here remains valid for polar molecules with an 
electronic spin in a strong magnetic field with the Zee- 
man splitting larger than the energy scales of the shield 
and the spin-rotation coupling; such a situation naturally 
appears for some polar molecules in magnetic traps. 

We choose to apply a static electric field E dc = E dc e z 
along the z-axis. For each molecule, a suitable basis set 
for the internal states is given by the eigenstates of the ro- 
tor Hamiltonian in the external static field. These states 
and the corresponding energies depend on the dimcn- 
sionlcss parameter dE dc / B and are denoted by \j, m z }i 
and Ej m „ , respectively, with m z the angular momen- 
tum along the z-axis and j denoting the different en- 
ergy manifolds, see Fig. QJa). I n addition, we apply a 
circularly (<r+) polarized microwave field E ac (t) prop- 
agating along the z-axis and coupling dominantly the 
ground state |0, 0)^ with the first excited state 
The microwave field is characterized by the detuning 
A = u) — — Eofi)/h and Rabi frequency ft = E^d c /h 
with the dipole coupling d c = | (0, 0|dj|l, 1) |. The lead- 
ing effect of the microwave field on the internal structure 
of each molecule is to mix the ground state |0,0) with 
the excited state |1, 1). We are interested in the regime 
with A, ft <C B/H and dE dc < 2B, where the rotating 
wave approximation is valid. In the rotating frame, these 
dressed states then take the form |+) = a\0, 0) + /3\1, 1} 
and |— ) = /3|0, 0)— a\l, 1) with the energy splitting AE = 
HVA 2 + 4ft 2 and a = -A/V A 2 + ft 2 , p = ft/V^Tft 2 , 
and A = (A + VA 2 +4ft 2 )/2. Throughout this let- 
ter, we are interested in a shield with a high barrier, 
which is optimized for parameters close to dE dc /B = 1, 
HA = 0.0155, and ft/A = 0.9258 (27| (see below). 

We next turn to the dipolar interaction and derive 
the dressed Born-Opcnhcimer potentials. Each polar 
molecule is prepared into the internal state by an 
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FIG. 2: Born-Oppenheimer potentials in the limit r » r B : (a) 
9 = and (b) 9 = ir/2. The effective potential V otf (r) (solid 
line) is repuslive for all angles 9. The dotted line denotes the 
antisymmetric level relevant during a three-body collision. 



adiabatic switching on of the microwave field. Conse- 
quently, the effective interaction potential V a{f (r) is deter- 
mined by the dressed Born-Oppenheimer potential adia- 
batically connected to the state |+)i|+) 2 , see Fig. [2] The 
competition between the rotational splitting B and the 
dipolc-dipole interaction provides a characteristic length 
scale r B = (d 2 /B) 1 ^ 3 . For large interparticle distances 
r » r B , the dipolc-dipole interaction is weak and does not 
couple different rotor levels. Consequently, the only rele- 
vant coupling appears due to the microwave field bewteen 
the manifolds with j = and j = 1. The relevant inter- 
nal states are then given by the three states \g)i = |0, 0),-, 
\e)i = |1, and \e)i = |1, —1);. As the microwave field 
couples each polar molecule with the same phase, the 
Born-Oppenheimer potentials separate into 6 symmetric 
and 3 antisymmetric dressed potentials. The effective 
potential adiabatically connected to the state |+)i|+)2 
is symmetric, and therefore, we can restrict the analysis 
to the symmeric potentials: a basis is given by the sym- 
metric states \g,g), \e,g), \e,e), \g,e), |e, e), and |e, e). 
Within the rotating frame, the Hamiltonian projected 
onto this subspace reduces to 
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with the dipole moments d g = \(g\di\g)\, d e = |(e|dj|e)|, 
and H eg = (d e d g - d 2 c /2)u - ^A and H ee = d\v - 2/iA. 
The terms v = (1 - 3 cos 2 d)/r 3 and fx = 3 sin 2 9e 2 ^/r 3 
describe the spatial dependence of the dipolc-dipole in- 
teraction, with 8 and <fi being the polar and azimuthal an- 
gles of r, respectively The Born-Oppenheimer potentials 
then follow from a diagonalization of the Hamiltonian H 
and are shown in Fig. [21 with the level adiabatically con- 
nected to the state |+)i|+) 2 (solid line) giving rise to the 
effective interaction V^ ff (r). 

The detuning A introduces a new length scale in the 
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problem ?'a = (d 2 /hA) 1 / 3 r B , e.g.. for LiCs with 
B/h ss 5.8GHz, r B w 9.2 nm and ta ~ 37.5nm at 
hA = 0.0155. At large interparticle distances r > r^, 
the Born-Oppenheimer potentials are well described by 
perturbation theory in the dipole-dipole interaction. The 
static electric field gives rise to a dipole moment d z = 
a 2 d g + p 2 d e along the z-axis, while the microwave field 
drives an additional dipole moment d± = \f2a[3d c ro- 
tating with frequency uj in the x-y plane. The com- 
bination of the two dipole forces provides the interac- 
tion Kff(r) = {d 2 z - d\/2) (l - 3cos 2 6») /r 3 . A proper 
choice of the two parameters E dc d/B and fi/A gives 
d z = d±/i/2, providing a cancelation of the leading 
dipole-dipole interaction The remaining interaction 
then follows from second order perturbation theory and 
provides a van-der-Waals-type repulsion 

Kb (r) = -4 [cf ' (1 - 3 cos 2 6) 2 + Cf ] 9 sin 4 o] (4) 



with 



hC^ 



^a 2 (3 2 [(d e -d g ) 2 +d 2 c ] 2 (5) 
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For the optimal values dE dc /B 



3VA 2 + 4ft 2 
1 and ft/A = 



0.9258, 



the van dcr Waals coefficients take the form C, 



0.004ftAr A and C 6 



(2) 
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< r A , 



(0) 
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0.005ftAr A . 



At shorter dis- 
thc effective interaction reduces 



tances r B <C r 
to V e!t (r) — (d 2 + d g d e [l — 3 cos 2 6\) jr 3 and remains re- 
pulsive for all angles 6. Thus it is possible to create purely 
repulsive interaction with large and adjustable strength. 

In order to determine the height of the potential bar- 
rier, a detailed analysis including all internal levels is 
required. Such a procedure is achieved by first deriv- 
ing Born-Oppenheimer potentials accounting for the cou- 
pling of the internal states \ j, m)i by the dipole-dipole in- 
teraction. In the second step, the microwave field, which 
couples these Born-Oppenheimer potentials, is included 
within a rotating wave approximation. The new dressed 
levels in the rotating frame are shown in Fig. [3Ja). It 
follows that the height of the shield (solid line) is limited 
by small avoided crossings. The first crossing (labeled 
A) appears with the level adiabatically connected to the 
symmetric state |1, 0; 0, 0) for a relative orientation of the 
molecules along the z-axis with ss 0, and it limits the 
barrier height of the shield to -E shio i d ps 0.02i? for the op- 
timal parameters. The radius R c for the breakdown of 
the shield is in the range R c ~ r B , which is still large 
compared to the distances where additional short range 
interactions have to be taken into account. This justifies 
our approach of restricting the analysis to the electric 
dipole-dipole interaction alone. 
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FIG. 3: (a) Born-Oppenheimer potentials for hA = 0.015B 
at 6 = 0. The plot shows all potentials (dashed) accessible via 
single-photon transitions from the state (solid) adiabatically 
connected to |+)i|+)2- The first crossing limiting the height 
of the shield appears in the region A for angles 6 ~ 0. (b) In- 
elastic cross section a ia due to diabatic crossings for different 
detunings A as a function of incoming kinetic energy Skin- 



Next, we analyze the validity of the Born-Oppenheimer 
approximation and study the influence of the kinetic en- 
ergy coupling different dressed potentials during a colli- 
sion. The influence of the kinetic energy is determined 
by the dimensionless parameter 7 = d 2 m//i 2 r B . For 
7 ^> 1, we can apply semiclassical theory; this condi- 
tion is well satisfied for typical polar molecules like LiCs 
with 7 w 6900. For a collision whose relative kinetic 
energy E kin is below the shield barrier, the processes 
giving rise to inelastic loss are (i) diabatic crossing be- 
tween different Born-Oppenheimer levels and (ii) quan- 
tum mechanical tunneling through the barrier. We start 
by studying the diabatic transitions first: the inelastic 
cross section within semi-classical approximation is com- 
puted by first determining the classical trajectory r cl (t) 
of a collision with a fixed impact parameter. To deter- 
mine the Landau-Zener diabatic crossings, we solve the 
full Schrodinger equation for the internal structure using 
the given relative motion r cl (t). The loss propability is 
then determined by the depiction of the adiabatic Born- 
Oppenheimer level at the classical turning point. Av- 
eraging over different impact parameters and angles of 
approach provides the inelastic cross section a ia due to 
diabatic transitions. The main contribution comes from 
the Born-Oppenheimer level closely approaching the ef- 
fective potential close to 6 w 0, see Figs. and[3Ja). 
Note that the standard Landau-Zener tunneling expres- 
sion can not be applied here as the levels have no real 
crossing. The inelastic cross section for different A and 
E Un is shown in Fig. 0(b) . We find an algebraic behavior 
of the inelastic cross section with cr in = p(E kill /hA) K r 2 
with re w 2.2 and p = 0.043 at hA = 0.015B, solid line 
in Fig. EIb). The loss rate l/r in during a two-particle 
collisions reduces to 
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where we have replaced the collision energy with the tem- 
perature T of the gas. Consequently, for LiCs at char- 
acteristic densities n ~ 10 13 cm -3 , the lifetime is several 
seconds even for T ~ lmK < HA. 

The second scenario for an inelastic collision is quan- 
tum mechanical tunneling through the barrier. For 7> 1 
such tunneling processes are strongly suppressed and can 
be studied using semiclassical techniques like WKB. The 
tunneling propability during a single collision is then 
given by the Eucledian action for the trajectory C start- 
ing at the classical turning point R$ and ending at the 
inner distance R c , where the shield starts to break down; 



dimensional crystalline phases. The many-body Hamil- 
tonian for a gas of ultracold polar molecules reduces to 



cxp 



ds yjm [V Bt{ (r) - E kiM ]/K 



(7) 



Thus, the highest tunneling probability is obtained along 
the collision axis with the lowest shield barrier, see discus- 
sion above. Then the characteritic scale for the tunneling 
amplitude at low incoming kinetic energies (i? kin <C hA) 
is given by P W kb = cxp(— c^/7), where the numerical con- 
stant for HA/B = 0.015 takes the form c sw 0.32. Conse- 
quently, for LiCs the tunneling is strongly suppressed and 
can be safely ignored at low kinetic energies E kin < HA. 

It is important to note that the qualitative behavior 
of the shield remains robust during a three-body colli- 
sion. The main modification to the shield is that the 
antisymmetric levels can open up an avoided crossing, as 
the parity symmetry present in the two-particle problem 
is broken for three particles. The relevant antisymmet- 
ric level is shown by a dotted line in Fig. [UJb). As this 
crossing appears at energies ~ftA, it does not modify the 
validity of the above discussion for incoming kinetic ener- 
gies E kin < HA. Thus, the shield prevents three particles 
from approaching each other on distances, where the for- 
mation of bound states can take place, and three-body 
losses are therefore effectively quenched. 

To summarize, we showed that properly adjusted con- 
tinous wave microwave and DC electric fields can create a 
repulsive shield resulting in large suppression of inelastic 
collisions. We now describe possible avenues opened by 
this work. For efficient evaporative cooling, it is impor- 
tant that elastic collisions allow for fast thcrmalization. 
The elastic scattering cross section at low collisional en- 
ergies is determined by the s-wave scattering length. In 
the present situation, the scattering length can be esti- 
mated via a partial wave expansion [25]: the dominant 
part is determined by the isotropic part of the effec- 
tive interaction potential vj® (r) = C^/r 6 . The s-wave 

1 /4 

scattering length then reduces to a s ~ (Cem/H 2 ) ~ 

ta (dPrn/h 2 r&) . For LiCs at the optimal detuning 
HA = 0.01 55, it follows that a s « 66 nm yielding a large 
elastic cross section ct o1 with suppressed losses providing 
an ideal system for evaporative cooling, e.g. a c] /a iD ~ 10 6 
for er in at T = 1 mK and er cl = 8ira 2 . 

Another application is the creation of three- 



c fi 



2 i r * 



(8) 



where we neglect additional terms due to the anisotropy 
of the effective potential and due to three-body inter- 
actions discussed in In analogy to the appearance 
of crystalline phases for polar molecules confined to two 
dimensions [261 ]. for strong repulsive van der Waals in- 
teractions, the system will undergo a phase transition 
from a liquid phase to a solid phase. The dimension- 
less parameter controlling the transition takes the form 



A = C 6 m/h 2 a 4 



ith 



1/a 3 the particle density: for 



weak interactions A <C 1 the ground state is in a liquid 
phase, while for strong interactions with A > 1 the sys- 
tem is characterized by a solid phase with broken transla- 
tional symmetry. Consequently, this opens up a way to- 
wards the creation of three-dimensional crystalline struc- 
tures with ultracold molecular gases. 
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